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I. INTRODUCTION 

If neutrinos turn out to be Majorana particles, their effects should manifest in lepton 
number violating (LNV) processes where the total lepton number changes in two units 
(AL = 2) [l, Q]. Conversely, the observation of decay/production phenomena with AL = 2 
would be very helpful in elucidating the mechanism of neutrino mass generation 3|. Given 
its relevance, it is very important to study all possible channels that may be sensitive to the 
effects of AL = 2 interactions and explore the constrains that they provide on the parameter 
space of specific models. 

In the case that LNV are induced by the exchange of Majorana neutrinos, their masses and 
mixing angles can be constrained from the experimental upper limits of AL = 2 observables. 
It is very well known that the most sensitive channels to very light Majorana neutrinos are 

. On the other hand, AL = 2 decays of 
ll| have proven to be useful to constrain 
sterile neutrinos with masses in the MeV to a few GeV range, which can be produced on 
their mass-s he,, a . the. .eacaon. Majo.a„a pa.acle. of th. Hnd are known as resonant 
neutrinos [5[. In addition to these widely studied three-body decays 



neutrinoless double beta decays of some nuclei J4 
pseudoscalar mesons jsl-ioi and tau leptons js], [gT 



m previous 



papers we have reported results on the analysis of four-body decays of neutral mesons [17 1 
and tau leptons 18| (see also 19|). It has been shown H, 20| that searches of these new 



(yet unexplored) decay channels can provide constraints on the parameter space of Majorana 
neutrinos that are complementary to three-body decays. 

In the present Brief Report we study the constraints that can be gotten from current 
experimental bounds on four-body decays of heavy mesons. Our study is motivated by 
searches reported recently by the LHCb collaboration, namely [l^ 



B{B- D\+fx-fx-) < 1.5 X 10" 



(1) 



This decay channel may receive contributions from the exchange of heavy Majorana neutri- 
nos via two different Feynman diagrams. We prove that in this case the Cabibbo-allowed 
channel, similar to the dominant contribution that underlies neutral B meson decays jl7|, 
is favored over the Cabibbo-suppressed one in most of the neutrino mass region. 

In addition, we also study the constrains provided by the four-body LNV decay channels 
of mesons. Althought experimental searches for these decays were reported long ago 
by the E791 col 



aboration 



2l| . no theoretical studies have been done so far. The reported 



upper limits 2l| currently are very mild, however our study shows that the constrains on 
the parameter space of Majorana neutrinos can become competitive and complementary to 
the ones gotten from three-body decays with improved limits on the branching fractions. 
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II. FOUR-BODY LNV DECAYS 



The Feynman diagrams that contribute to the LNV decay — t- D^tt^ fi^ fi^ are shown 
in Fig. [H As in previous studies, we assume that only one heavy neutrino A^, with a mass 
such that it can be produced resonantly in B decays, dominates the decay amp 



n 



itude. For 



such heavy neutrinos, the diagram of Figure Ic gives a negligible contribution [7|. 

The amplitude in Figure la can be resonantly enhanced for neutrino masses in the range 
ttItt + 171^ < rriN < iJiB — ttld — m^. With obvious notation for particle momenta, it can be 
written as 0, Q 



GlV^''''V^r{D{pD)\crb\B{pBm.pl]L^p^ 



(2) 




FIG. 1. Lowest order (a) spectator, (6) annihilation, and (c) t-channel diagrams, mediated by a 
heavy Majorana neutrino N in — )■ D^tt^ fi^ ^i" decay. 



where \V2™\ = 40.9x10^^ and \V^™\ = 0.97425 [12j are the Cabibbo-Kobayashi-Maskawa 



quark mixing matrix elements, Gp is the Fermi constant and = 130.4 MeV [1^ the pion 
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decay constant. The information about the Majorana neutrino exchange is included in the 
leptonic tensor Lafs (properly antisymmetrized under exchange of lepton momenta ^1^2) jsl 



Lal3 =\V^,N\'^mN U{pi) 



PBU''ip2), (3) 



with Pr = (1 + 75)72, Q = pi + P2 + Ptt the momentum transfer, rriN and Tn denote the 
mass and decay width of the heavy neutrino, respectively (F^v can be calculated as shown 
in Ref. jsl). The constant V^at denotes the mixing of the muon with the heavy neutrino N 
in the charged current interaction. 

The hadronic matrix element in Eq. (|2]) is given by 



{D{pn)\c^ab\B{pB)) =[{pB+PD)a - J Qo)F^^''{t) + -g„Fo^-^(t), 



(4) 



where A = (m^ — Fq^^ are the scalar and vector form factors for the B ^ D 

transition evaluated aX t = Q"^ . We will use the predictions for these form factors that are 
obtained from Lattice QCD calculations 22 1. 

The contribution from Figure lb is similar to the one encountered in three-body decays 
— 7- X^fi~fi~ with the subsequent decay of the charmed resonance — t- D^Tr~^. The 
corresponding amplitude becomes resonantly enhanced for neutrino mass values in the range 
rriu + + < m^r < itlb — m^. The associated decay amplitude is given by 



(5) 



where IK™ I = 4.15 x IQ-^, |V;^™| = 0.230 [jj, /b = 194 MeV [jj and the leptonic 
tensor 



Lal3 = \Vf,N\'^mN u{pi) 



7a 7/3 



(Pb - Pi 



)2 - m% + iVNTriN 



+ (a O 13, pi 



Pru'{p2)- (6) 



The hadronic matrix element is parametrized as follows 

A'. 



{D{pD)T^{p^)\h^c\Q) = [{pD-PnT + ^k-)F['~'^{e) + ^k'^F.^^^ik 



(7) 



^. The form factors will be modeled as F['f^'^{k'^ 



where k = pn + p-j^ and A' = — m 
F{^^^iO)-BWxSk^), with FiV^(O) = 0.67 [23] its value at zero momentum transfer. We will 
use a simplified model where the Breit-Wigner (BW) function introduced above is dominated 
by a single resonance: 



BWx^k' 



m 



X, 



I He. 



k"^ - imxTx, 



(8) 



where the mass and width of the charmed resonances X^. correspond to the Z)*"'"(2010) and 
1^0^(2400), respectively, for the spin-1 and spin-0 form factors. 
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Althought the annihilation amphtude in Figure lb is Cabibbo-suppressed with respect 
to the spectator amphtude from Fig. la, a strong enhancement of the former due to the 
exchange of the narrow D*(2010) resonance is possible. Note that the sensitivity to the mass 
of Majorana neutrinos explored by diagrams la and lb are complementary and overlap only 
over a small window 2.1 GeV < rriN < 3.3 GeV. 



The phase space of the 



cinematical variables 



24 



bur-body decay is determined by the limits on the five independent 



25| . We perform the numerical integrations using the VEGAS code 



261] and implementing the single-diagram enhanced channel integration method 27|]. Figure 



[2]shows the excluded region (above the solid and short-dashed curves) for | V^at^ as a function 
of m^v that is obtained from the upper limit reported by the LHCb Collaboration, Eq. ([1]). 
As it was pointed out above, the dominant effects of Majorana neutrinos come the diagram 
in Figure la while the effects from diagram lb are visible only at higher values of rriN- For 
comparison, we also display the updated exclusion plots obtained from searches of three-body 
{D^ , , B^) — !■ 7r+/i^yU^ decays. For neutrino mass values below 1 GeV, the four-body 
decay under consideration is able to exclude a larger region of the |V^Arp than the Cabibbo- 
suppressed channel — )• n^fi^fi^ (long-dashed plot in Figure 2) despite the fact that the 



upper limit on the later is of order 10 ^ IJ]. Therefore, future and improved upper limits on 



the B{B — i- D^Ti^fi /i ) will be very useful to derive better constrains on the fiN mixing 
angle. 




FIG. 2. Exclusion regions on heavy mixing |V^Arp as a function of the Majorana neutrino mass 
rriN, from LNV dimuon modes of heavy mesons. The solid and short-dashed 

An analogous process to the one studied before is the — t- Tr'^vr^/i^yU^ decay. If this 
AL = 2 decay is due to the exchange of Majorana neutrinos, we should have similar contri- 
butions to the amplitude as the ones shown in Figure 1. The corresponding spectator and an- 
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nihilation contributions of Figs la and lb are both proportional to the product V^^^V^^, 
therefore we can expect that the p~^{770) resonance intermediate contribution (Fig. lb) 
becomes important in this case. Moreover, the sensitivity region of neutrino masses where 
spectator and annihilation amplitudes are resonantly enhanced are very similar. However, 
since this four-body decay has similar Cabibbo-suppressed couplings as — t- 7i^fi^fi~ we 
expect that the excluded region on Majorana neutrinos from experimental upper limits on 
both decays would be very similar. 



III. FOUR-BODY LNV DECAYS 



The four-body D° — )■ h^h'^l^l'^ decays, where h, h' = n, K and /, /' = e, fi, were studied 
by the E791 collaboration [21i] more than a decade ago and 90% CL upper limits of order 
10~^ ~ 10^^ for their branching ratios were obtained. In this section we focus on the 

—7- {K^Ti^ , 7i^7i^)fi^ channels in order to illustrate the potential these decays offer to 
constrain the parameters of Majorana neutrinos. The upper limits for these decays obtained 
in Ref. 2l| are: 

B{D° n-n-fi+fi^) < 2.9 x 10"^ 

^ K-n-fi+i2+) < 3.9 X 10"^ (9) 

Since weak currents carry one unit of electric charge, the Feynman diagram analogous 
to Figure lb does not contribute in this case. The dominant contribution to the decay 
amplitude is given by the spectator diagram similar to Figure la and its amplitude becomes: 

-MS = G'^Kf ''KF™(MP/.)|g7"c|D(p^))[z/.p^]L„,, (10) 

with q = d{s) for h = -n{K). The leptonic tensor current Lap is given by Eq. ([3D and the 

For the purposes 
In the case 



hadronic current (/i(ph)|^7°c|-D(pD)) can be parametrized as in Eq. ( 
of our numerical evaluation we take the form factors F^^^ from Ref. 
of identical mesons in the final state the decay amplitude must, in addition, be symmetrized 
under exchange of their momenta. 

By using the upper limits given in ([9]), in Fig. [3] we plot the exclusion regions provided by 
—7- TT^TT^/i^/i^ (solid line) and — > K^tt^ p^ p^ (short-dashed line) decays. Currently, 
these decays provide only loose constrains on Majorana neutrinos given the poor upper 
limits available on branching ratios. Note however that improvements by two or three 
orders of magnitude on the branching ratios would yield competetive constraints on the 
mixing angle |V"^Arp compared to other LNV meson decays, as it can be appreciated in 
Figure [31 This goal is certainly at the reach of current (LHCb, BESIII) and future (Belle 
II) D meson factories. 



10 = 



10' 



2x10"' 3x10"' 



3 4 5 
[GeV] 



FIG. 3. Exclusion curves for |V^Arp as a function of mjv from four-body LNV decays (solid 
and short-dashed lines). Constraints from other heavy meson decays are shown for comparison. 



Summarizing, in this Report we derive the first constraints on parameters of Majo- 
rana neutrinos that can be obtained from the LNV decays — t- D^tt^ and — >■ 
(i^"7r^/7r^7r^)/i+/i+. We use the upper hmits on the corresponding branching ratios that 
have been reported, respectively, by the LHCb [l^ and the E791 [21| Collaborations. As in 
previous works, we assume that only one Majorana neutrino A^, with a mass of hundreds of 
MeV to a few GeV, contribute resonantly to the decay amphtude (LNV heavy meson decays 
are insensitive to either, very light or very heavy neutrinos 29|). The excluded region in the 
iV^ATp vs. niN plane obtained from the upper limit on B{B~ — )• D^-n^ jj,^ ji^) reported by 
LHCb [l^ is similar to the one obtained from existing bounds on LNV three body decays 
of charged D and B mesons. This is possible because the four-body decay is a Cabibbo- 
allowed decay, compared to the three-body decay which proceeds via a Cabibbo-suppressed 
mechanism. 



We have also studied the four-body LNV decays of the neutral charmed meson. We use 
the two most restrictive upper bounds obtained by the E791 collaboration to derive bounds 
on the same parameters of Majorana neutrinos. These bounds turn out to be very mild 
at present, but their improvements at current and planned D meson factories would make 
them competitive with constrains derived from three-body decays of charged mesons. 
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